Spatial Architecture of Myeloid and T Cells Orchestrates Immune Evasion and Clinical Outcome in Lung Cancer.

Understanding the role of the tumour microenvironment (TME) in lung cancer is critical to improving patient outcome. We identified four histology-independent archetype TMEs in treatment-naive early-stage lung cancer using imaging mass cytometry in the TRACERx study (n=81 patients/198 samples/2.3million cells). In immune-hot adenocarcinomas, spatial niches of T cells and macrophages increased with clonal neoantigen burden, whereas such an increase was observed for niches of plasma and B cells in immune-excluded squamous cell carcinomas (LUSC). Immune-low TMEs were associated with fibroblast barriers to immune infiltration. The fourth archetype, characterised by sparse lymphocytes and high tumour-associated neutrophil (TAN) infiltration, had tumour cells spatially separated from vasculature and exhibited low spatial intratumour heterogeneity. TAN-High LUSC had frequent PIK3CA mutations. TAN-High tumours harboured recently expanded and metastasis-seeding subclones and had a shorter disease-free survival independent of stage. These findings delineate genomic, immune and physical barriers to immune surveillance and implicate neutrophil-rich TMEs in metastasis.

Human endogenous retrovirus onco-exaptation counters cancer cell senescence through calbindin.

Increased levels and diversity of human endogenous retrovirus (HERV) transcription characterize most cancer types and are linked with disease outcomes. However, the underlying processes are incompletely understood. Here, we show that elevated transcription of HERVH proviruses predicted survival of lung squamous cell carcinoma (LUSC) and identified an isoform of CALB1, encoding calbindin, ectopically driven by an upstream HERVH provirus under the control of KLF5, as the mediator of this effect. HERVH-CALB1 expression was initiated in preinvasive lesions and associated with their progression. Calbindin loss in LUSC cell lines impaired in vitro and in vivo growth and triggered senescence, consistent with a protumor effect. However, calbindin also directly controlled the senescence-associated secretory phenotype (SASP), marked by secretion of CXCL8 and other neutrophil chemoattractants. In established carcinomas, CALB1-negative cancer cells became the dominant source of CXCL8, correlating with neutrophil infiltration and worse prognosis. Thus, HERVH-CALB1 expression in LUSC may display antagonistic pleiotropy, whereby the benefits of escaping senescence early during cancer initiation and clonal competition were offset by the prevention of SASP and protumor inflammation at later stages.

Evolutionary characterization of lung adenocarcinoma morphology in TRACERx.

Lung adenocarcinomas (LUADs) display a broad histological spectrum from low-grade lepidic tumors through to mid-grade acinar and papillary and high-grade solid, cribriform and micropapillary tumors. How morphology reflects tumor evolution and disease progression is poorly understood. Whole-exome sequencing data generated from 805 primary tumor regions and 121 paired metastatic samples across 248 LUADs from the TRACERx 421 cohort, together with RNA-sequencing data from 463 primary tumor regions, were integrated with detailed whole-tumor and regional histopathological analysis. Tumors with predominantly high-grade patterns showed increased chromosomal complexity, with higher burden of loss of heterozygosity and subclonal somatic copy number alterations. Individual regions in predominantly high-grade pattern tumors exhibited higher proliferation and lower clonal diversity, potentially reflecting large recent subclonal expansions. Co-occurrence of truncal loss of chromosomes 3p and 3q was enriched in predominantly low-/mid-grade tumors, while purely undifferentiated solid-pattern tumors had a higher frequency of truncal arm or focal 3q gains and SMARCA4 gene alterations compared with mixed-pattern tumors with a solid component, suggesting distinct evolutionary trajectories. Clonal evolution analysis revealed that tumors tend to evolve toward higher-grade patterns. The presence of micropapillary pattern and 'tumor spread through air spaces' were associated with intrathoracic recurrence, in contrast to the presence of solid/cribriform patterns, necrosis and preoperative circulating tumor DNA detection, which were associated with extra-thoracic recurrence. These data provide insights into the relationship between LUAD morphology, the underlying evolutionary genomic landscape, and clinical and anatomical relapse risk.

Antibodies against endogenous retroviruses promote lung cancer immunotherapy.

B cells are frequently found in the margins of solid tumours as organized follicles in ectopic lymphoid organs called tertiary lymphoid structures (TLS)1,2. Although TLS have been found to correlate with improved patient survival and response to immune checkpoint blockade (ICB), the underlying mechanisms of this association remain elusive1,2. Here we investigate lung-resident B cell responses in patients from the TRACERx 421 (Tracking Non-Small-Cell Lung Cancer Evolution Through Therapy) and other lung cancer cohorts, and in a recently established immunogenic mouse model for lung adenocarcinoma3. We find that both human and mouse lung adenocarcinomas elicit local germinal centre responses and tumour-binding antibodies, and further identify endogenous retrovirus (ERV) envelope glycoproteins as a dominant anti-tumour antibody target. ERV-targeting B cell responses are amplified by ICB in both humans and mice, and by targeted inhibition of KRAS(G12C) in the mouse model. ERV-reactive antibodies exert anti-tumour activity that extends survival in the mouse model, and ERV expression predicts the outcome of ICB in human lung adenocarcinoma. Finally, we find that effective immunotherapy in the mouse model requires CXCL13-dependent TLS formation. Conversely, therapeutic CXCL13 treatment potentiates anti-tumour immunity and synergizes with ICB. Our findings provide a possible mechanistic basis for the association of TLS with immunotherapy response.

Spatial Positioning of Immune Hotspots Reflects the Interplay between B and T Cells in Lung Squamous Cell Carcinoma.

Beyond tertiary lymphoid structures, a significant number of immune-rich areas without germinal center-like structures are observed in non-small cell lung cancer. Here, we integrated transcriptomic data and digital pathology images to study the prognostic implications, spatial locations, and constitution of immune rich areas (immune hotspots) in a cohort of 935 patients with lung cancer from The Cancer Genome Atlas. A high intratumoral immune hotspot score, which measures the proportion of immune hotspots interfacing with tumor islands, was correlated with poor overall survival in lung squamous cell carcinoma but not in lung adenocarcinoma. Lung squamous cell carcinomas with high intratumoral immune hotspot scores were characterized by consistent upregulation of B-cell signatures. Spatial statistical analyses conducted on serial multiplex IHC slides further revealed that only 4.87% of peritumoral immune hotspots and 0.26% of intratumoral immune hotspots were tertiary lymphoid structures. Significantly lower densities of CD20+CXCR5+ and CD79b+ B cells and less diverse immune cell interactions were found in intratumoral immune hotspots compared with peritumoral immune hotspots. Furthermore, there was a negative correlation between the percentages of CD8+ T cells and T regulatory cells in intratumoral but not in peritumoral immune hotspots, with tertiary lymphoid structures excluded. These findings suggest that the intratumoral immune hotspots reflect an immunosuppressive niche compared with peritumoral immune hotspots, independent of the distribution of tertiary lymphoid structures. A balance toward increased intratumoral immune hotspots is indicative of a compromised antitumor immune response and poor outcome in lung squamous cell carcinoma. SIGNIFICANCE: Intratumoral immune hotspots beyond tertiary lymphoid structures reflect an immunosuppressive microenvironment, different from peritumoral immune hotspots, warranting further study in the context of immunotherapies.

Characterisation of tumour microenvironment remodelling following oncogene inhibition in preclinical studies with imaging mass cytometry.

Mouse models are critical in pre-clinical studies of cancer therapy, allowing dissection of mechanisms through chemical and genetic manipulations that are not feasible in the clinical setting. In studies of the tumour microenvironment (TME), multiplexed imaging methods can provide a rich source of information. However, the application of such technologies in mouse tissues is still in its infancy. Here we present a workflow for studying the TME using imaging mass cytometry with a panel of 27 antibodies on frozen mouse tissues. We optimise and validate image segmentation strategies and automate the process in a Nextflow-based pipeline (imcyto) that is scalable and portable, allowing for parallelised segmentation of large multi-image datasets. With these methods we interrogate the remodelling of the TME induced by a KRAS G12C inhibitor in an immune competent mouse orthotopic lung cancer model, highlighting the infiltration and activation of antigen presenting cells and effector cells.

Assessment of long non-coding RNA expression reveals novel mediators of the lung tumour immune response.

The tumour immune microenvironment is a crucial mediator of lung tumourigenesis, and characterizing the immune landscape of patient tumours may guide immunotherapy treatment regimens and uncover novel intervention points. We sought to identify the landscape of tumour-infiltrating immune cells in the context of long non-coding RNA (lncRNAs), known regulators of gene expression. We examined the lncRNA profiles of lung adenocarcinoma (LUAD) tumours by interrogating RNA sequencing data from microdissected and non-microdissected samples (BCCRC and TCGA). Subsequently, analysis of single-cell RNA sequencing data from lung tumours and flow-sorted healthy peripheral blood mononuclear cells identified lncRNAs in immune cells, highlighting their biological and prognostic relevance. We discovered lncRNA expression patterns indicative of regulatory relationships with immune-related protein-coding genes, including the relationship between AC008750.1 and NKG7 in NK cells. Activation of NK cells in vitro was sufficient to induce AC008750.1 expression. Finally, siRNA-mediated knockdown of AC008750.1 significantly impaired both the expression of NKG7 and the anti-tumour capacity of NK cells. We present an atlas of cancer-cell extrinsic immune cell-expressed lncRNAs, in vitro evidence for a functional role of lncRNAs in anti-tumour immune activity, which upon further exploration may reveal novel clinical utility as markers of immune infiltration.

Deregulation of a Cis-Acting lncRNA in Non-small Cell Lung Cancer May Control HMGA1 Expression.

BACKGROUND: Long non-coding RNAs (lncRNAs) have long been implicated in cancer-associated phenotypes. Recently, a class of lncRNAs, known as cis-acting, have been shown to regulate the expression of neighboring protein-coding genes and may represent undiscovered therapeutic action points. The chromatin architecture modification gene HMGA1 has recently been described to be aberrantly expressed in lung adenocarcinoma (LUAD). However, the mechanisms mediating the expression of HMGA1 in LUAD remain unknown. Here we investigate the deregulation of a putative cis-acting lncRNA in LUAD, and its effect on the oncogene HMGA1. METHODS: LncRNA expression was determined from RNA-sequencing data of tumor and matched non-malignant tissues from 36 LUAD patients. Transcripts with significantly deregulated expression were identified and validated in a secondary LUAD RNA-seq dataset (TCGA). SiRNA-mediated knockdown of a candidate cis-acting lncRNA was performed in BEAS-2B cells. Quantitative real-time PCR was used to observe the effects of lncRNA knockdown on the expression of HMGA1. RESULTS: We identified the lncRNA RP11.513I15.6, which we refer to as HMGA1-lnc, neighboring HMGA1 to be significantly downregulated in both LUAD cohorts. Conversely, we found HMGA1 significantly overexpressed in LUAD and anticorrelated with HMGA1-lnc. In vitro experiments demonstrated siRNA-mediated inhibition of HMGA1-lnc in immortalized non-malignant lung epithelial cells resulted in a significant increase in HMGA1 gene expression. CONCLUSION: Our results suggest that HMGA1-lnc is a novel cis-acting lncRNA that negatively regulates HMGA1 gene expression in lung cells. Further characterization of this regulatory mechanism may advance our understanding of the maintenance of lung cancer phenotypes and uncover a novel therapeutic intervention point for tumors driven by HMGA1.

Epithelial tumor suppressor ELF3 is a lineage-specific amplified oncogene in lung adenocarcinoma.

Gene function in cancer is often cell type-specific. The epithelial cell-specific transcription factor ELF3 is a documented tumor suppressor in many epithelial tumors yet displays oncogenic properties in others. Here, we show that ELF3 is an oncogene in the adenocarcinoma subtype of lung cancer (LUAD), providing genetic, functional, and clinical evidence of subtype specificity. We discover a region of focal amplification at chromosome 1q32.1 encompassing the ELF3 locus in LUAD which is absent in the squamous subtype. Gene dosage and promoter hypomethylation affect the locus in up to 80% of LUAD analyzed. ELF3 expression was required for tumor growth and a pan-cancer expression network analysis supports its subtype and tissue specificity. We further show that ELF3 displays strong prognostic value in LUAD but not LUSC. We conclude that, contrary to many other tumors of epithelial origin, ELF3 is an oncogene and putative therapeutic target in LUAD.

Aberrant Expression of Pseudogene-Derived lncRNAs as an Alternative Mechanism of Cancer Gene Regulation in Lung Adenocarcinoma.

Transcriptome sequencing has led to the widespread identification of long non-coding RNAs (lncRNAs). Subsequently, these genes have been shown to hold functional importance in human cellular biology, which can be exploited by tumors to drive the hallmarks of cancer. Due to the complex tertiary structure and unknown binding motifs of lncRNAs, there is a growing disparity between the number of lncRNAs identified and those that have been functionally characterized. As such, lncRNAs deregulated in cancer may represent critical components of cancer pathways that could serve as novel therapeutic intervention points. Pseudogenes are non-coding DNA sequences that are defunct relatives of their protein-coding parent genes but retain high sequence similarity. Interestingly, certain lncRNAs expressed from pseudogene loci have been shown to regulate the protein-coding parent genes of these pseudogenes in trans particularly because of this sequence complementarity. We hypothesize that this phenomenon occurs more broadly than previously realized, and that aberrant expression of lncRNAs overlapping pseudogene loci provides an alternative mechanism of cancer gene deregulation. Using RNA-sequencing data from two cohorts of lung adenocarcinoma, each paired with patient-matched non-malignant lung samples, we discovered 104 deregulated pseudogene-derived lncRNAs. Remarkably, many of these deregulated lncRNAs (i) were expressed from the loci of pseudogenes related to known cancer genes, (ii) had expression that significantly correlated with protein-coding parent gene expression, and (iii) had lncRNA protein-coding parent gene expression that was significantly associated with survival. Here, we uncover evidence to suggest the lncRNA-pseudogene-protein-coding gene axis as a prominent mechanism of cancer gene regulation in lung adenocarcinoma, and highlights the clinical utility of exploring the non-coding regions of the cancer transcriptome.

Hyperspectral cell sociology reveals spatial tumor-immune cell interactions associated with lung cancer recurrence.

BACKGROUND: The tumor microenvironment (TME) is a complex mixture of tumor epithelium, stroma and immune cells, and the immune component of the TME is highly prognostic for tumor progression and patient outcome. In lung cancer, anti-PD-1 therapy significantly improves patient survival through activation of T cell cytotoxicity against tumor cells. Direct contact between CD8+ T cells and target cells is necessary for CD8+ T cell activity, indicating that spatial organization of immune cells within the TME reflects a critical process in anti-tumor immunity. Current immunohistochemistry (IHC) imaging techniques identify immune cell numbers and densities, but lack assessment of cell-cell spatial relationships (or "cell sociology"). Immune functionality, however, is often dictated by cell-to-cell contact and cannot be resolved by simple metrics of cell density (for example, number of cells per mm2). To address this issue, we developed a Hyperspectral Cell Sociology technology platform for the analysis of cell-cell interactions in multi-channel IHC-stained tissue. METHODS: Tissue sections of primary tumors from lung adenocarcinoma patients with known clinical outcome were stained using multiplex IHC for CD3, CD8, and CD79a, and hyperspectral image analysis determined the phenotype of all cells. A Voronoi diagram for each cell was used to approximate cell boundaries, and the cell type of all neighboring cells was identified and quantified. Monte Carlo analysis was used to assess whether cell sociology patterns were likely due to random distributions of the cells. RESULTS: High density of intra-tumoral CD8+ T cells was significantly associated with non-recurrence of tumors. A cell sociology pattern of CD8+ T cells surrounded by tumor cells was more significantly associated with non-recurrence compared to CD8+ T cell density alone. CD3+ CD8- T cells surrounded by tumor cells was also associated with non-recurrence, but at a similar significance as cell density alone. Cell sociology metrics improved recurrence classifications of 12 patients. Monte Carlo re-sampling analysis determined that these cell sociology patterns were non-random. CONCLUSION: Hyperspectral Cell Sociology expands our understanding of the complex interplay between tumor cells and immune infiltrate. This technology could improve predictions of responses to immunotherapy and lead to a deeper understanding of anti-tumor immunity.

Somatic mutation-associated T follicular helper cell elevation in lung adenocarcinoma.

T follicular helper cells (Tfh) play crucial roles in the development of humoral immunity. In the B cell-rich germinal center of lymphoid organs, they select for high-affinity B cells and aid in their maturation. While Tfh have known roles in B cell malignancies and have prognostic value in some epithelial cancers, their role in lung tumour initiation and development is unknown. Through immune cell deconvolution, we observed significantly increased Tfh in tumours from two independent cohorts of lung adenocarcinomas and found that this upregulation occurs early in tumour development. A subset of tumours were stained for T and B cells using multicolour immunohistochemistry, which revealed the presence of tumour-adjacent tertiary lymphoid organs in 17/20 cases each with an average of 16 Tfh observed in the germinal center. Importantly, Tfh levels were correlated with tumour mutational load and immunogenic cancer testis antigens, suggesting their involvement in mounting an active immune response against tumour neoantigens.

Expanding the miRNA Transcriptome of Human Kidney and Renal Cell Carcinoma.

Despite advancements in therapeutic strategies, diagnostic and prognostic molecular markers of kidney cancer remain scarce, particularly in patients who do not harbour well-defined driver mutations. Recent evidence suggests that a large proportion of the human noncoding transcriptome has escaped detection in early genomic explorations. Here, we undertake a large-scale analysis of small RNA-sequencing data from both clear cell renal cell carcinoma (ccRCC) and nonmalignant samples to generate a robust set of miRNAs that remain unannotated in kidney tissues. We find that these novel kidney miRNAs are also expressed in renal cancer cell lines. Moreover, these sequences are differentially expressed between ccRCC and matched nonmalignant tissues, implicating their involvement in ccRCC biology and potential utility as tumour-specific markers of disease. Indeed, we find some of these miRNAs to be significantly associated with patient survival. Finally, target prediction and subsequent pathway analysis reveals that miRNAs previously unannotated in kidney tissues may target genes involved in ccRCC tumourigenesis and disease biology. Taken together, our results represent a new resource for the study of kidney cancer and underscore the need to characterize the unexplored areas of the transcriptome.

Large-scale discovery of previously undetected microRNAs specific to human liver.

MicroRNAs (miRNAs) are crucial regulators of gene expression in normal development and cellular homeostasis. While miRNA repositories contain thousands of unique sequences, they primarily contain molecules that are conserved across several tissues, largely excluding lineage and tissue-specific miRNAs. By analyzing small non-coding RNA sequencing data for abundance and secondary RNA structure, we discovered 103 miRNA candidates previously undescribed in liver tissue. While expression of some of these unannotated sequences is restricted to non-malignant tissue, downregulation of most of the sequences was detected in liver tumors, indicating their importance in the maintenance of liver homeostasis. Furthermore, target prediction revealed the involvement of the unannotated miRNA candidates in fatty-acid metabolism and tissue regeneration, which are key pathways in liver biology. Here, we provide a comprehensive analysis of the undiscovered liver miRNA transcriptome, providing new resources for a deeper exploration of organ-specific biology and disease.

In Vivo Validation of PAPSS1 (3'-phosphoadenosine 5'-phosphosulfate synthase 1) as a Cisplatin-sensitizing Therapeutic Target.

Purpose: Our previous screening efforts found that inhibition of PAPSS1 increases the potency of DNA-damaging agents in non-small cell lung cancer (NSCLC) cell lines. Here, we explored the clinical relevance of PAPSS1 and further investigated it as a therapeutic target in preclinical model systems.Experimental Design: PAPSS1 expression and cisplatin IC50 values were assessed in 52 lung adenocarcinoma cell lines. Effects of PAPSS1 inhibition on A549 cisplatin sensitivity under hypoxic and starvation conditions, in 3D spheroids, as well as in zebrafish and mouse xenografts, were evaluated. Finally, the association between PAPSS1 expression levels and survival in patients treated with standard chemotherapy was assessed.Results: Our results show a positive correlation between low PAPSS1 expression and increased cisplatin sensitivity in lung adenocarcinoma. In vitro, the potentiation effect was greatest when A549 cells were serum-starved under hypoxic conditions. When treated with low-dose cisplatin, PAPSS1-deficient A549 spheroids showed a 58% reduction in size compared with control cells. In vivo, PAPSS1 suppression and low-dose cisplatin treatment inhibited proliferation of lung tumor cells in zebrafish xenografts and significantly delayed development of subcutaneous tumors in mice. Clinical data suggest that NSCLC and ovarian cancer patients with low PAPSS1 expression survive longer following platinum-based chemotherapy.Conclusions: These results suggest that PAPSS1 inhibition enhances cisplatin activity in multiple preclinical model systems and that low PAPSS1 expression may serve as a biomarker for platin sensitivity in cancer patients. Developing strategies to target PAPSS1 activity in conjunction with platinum-based chemotherapy may offer an approach to improving treatment outcomes. Clin Cancer Res; 23(21); 6555-66. ©2017 AACR.

DLK1-DIO3 imprinted locus deregulation in development, respiratory disease, and cancer.

INTRODUCTION: The imprinted DLK1-DIO3 locus at 14q32.1-32.31 holds biological significance in fetal development, whereby imprinting errors are causal to developmental disorders. Emerging evidence has implicated this locus in other diseases including cancer, highlighting the biological parallels between fetal organ and tumour development. Areas covered: Controlled regulation of gene expression from the imprinted DLK1-DIO3 locus at 14q32.1-32.31 is crucial for proper fetal development. Deregulation of locus gene expression due to imprinting errors has been mechanistically linked to the developmental disorders Kagami-Ogata Syndrome and Temple Syndrome. In adult tissues, deregulation of locus genes has been associated with multiple malignancies although the causal genetic mechanisms remain largely uncharacterised. Here, we summarize the genetic mechanisms underlying the developmental disorders that arise as a result of improper locus imprinting and the resulting developmental phenotypes, emphasizing both the coding and noncoding components of the locus. We further highlight biological parallels common to both fetal development and disease, with a specific focus on lung development, respiratory disease, and lung cancer. Expert commentary: Many commonalities between respiratory and developmental defects have emerged with respect to the 14q32 locus, emphasizing the importance of studying the effects of imprinting on gene regulation patterns at this locus in both biological settings.

Deregulation of small non-coding RNAs at the DLK1-DIO3 imprinted locus predicts lung cancer patient outcome.

Deregulation of the imprinted DLK1-DIO3 locus at chromosome 14q32.1-14q32.31 has been associated with developmental and respiratory disorders, including cancer. In lung cancer, deregulation of imprinting at DLK1-DIO3 was recently described in smokers. Deregulated expression of a microRNA (miRNA) cluster mapping to this locus was also associated with patient outcome, suggesting the importance of this locus to lung cancer disease phenotypes. The DLK1-DIO3 locus is complex, and encodes several protein-coding genes, in addition to long and short non-coding RNAs. While the role of miRNAs is established, the biological importance of another relevant class of small RNAs, PIWI-interacting RNAs (piRNAs), has not been investigated. When somatically expressed, piRNAs regulate gene transcription through DNA methylation. Interestingly, their expression patterns have been observed to be altered in cancer and correlated with patient outcome. Here, we characterize the somatic expression of piRNAs encoded at DLK1-DIO3 in two independent cohorts of lung adenocarcinoma and lung squamous cell carcinoma and investigate their associations with patient outcome. We find that the expression of piRNAs encoded at DLK1-DIO3 enhances the prognostic potential of small non-coding RNAs specific to this locus in predicting patient outcome, further emphasizing the importance of regulation at this locus in lung cancer.

Emerging roles of T helper 17 and regulatory T cells in lung cancer progression and metastasis.

Lung cancer is a leading cause of cancer-related deaths worldwide. Lung cancer risk factors, including smoking and exposure to environmental carcinogens, have been linked to chronic inflammation. An integral feature of inflammation is the activation, expansion and infiltration of diverse immune cell types, including CD4+ T cells. Within this T cell subset are immunosuppressive regulatory T (Treg) cells and pro-inflammatory T helper 17 (Th17) cells that act in a fine balance to regulate appropriate adaptive immune responses.In the context of lung cancer, evidence suggests that Tregs promote metastasis and metastatic tumor foci development. Additionally, Th17 cells have been shown to be an integral component of the inflammatory milieu in the tumor microenvironment, and potentially involved in promoting distinct lung tumor phenotypes. Studies have shown that the composition of Tregs and Th17 cells are altered in the tumor microenvironment, and that these two CD4+ T cell subsets play active roles in promoting lung cancer progression and metastasis.We review current knowledge on the influence of Treg and Th17 cells on lung cancer tumorigenesis, progression, metastasis and prognosis. Furthermore, we discuss the potential biological and clinical implications of the balance among Treg/Th17 cells in the context of the lung tumor microenvironment and highlight the potential prognostic function and relationship to metastasis in lung cancer.

HPV status is associated with altered PIWI-interacting RNA expression pattern in head and neck cancer.

OBJECTIVES: As HPV-induced cases of oral malignancy increase, it is important to understand the molecular differences between HPV positive and negative head and neck squamous cell carcinoma (HNSCC). PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs aberrantly expressed in cancer. We analyzed global piRNA expression patterns to define the HNSCC piRNA transcriptome and assess whether HPV infection status associates with changes in piRNA levels. MATERIALS AND METHODS: A total of 498 HNSCC small RNA sequencing libraries were acquired from the Cancer Genomics Hub (cgHUB) Data Repository and a custom sequence analysis pipeline was developed to deduce piRNA expression from raw sequencing data. Expression matrices were aligned to clinicopathological features in order to analyze piRNA expression patterns across different HNSCC groups. The association of a piRNA signature with HPV-positive patient survival was evaluated using a Cox proportional hazard model. RESULTS: Analysis of piRNA levels between HNSCC and non-malignant tissues revealed distinct expression patterns, with 87 piRNAs exclusively expressed in tumor samples. HPV infection status affected the expression of 41 of these piRNAs. Eleven (26.8%) piRNAs were significantly downregulated in HPV16/18 tumors compared to other HPV types. Remarkably, expression of a combination of five-piRNAs in HPV-positive HNSCC tumors was associated with worse overall survival. CONCLUSION: The expression of specific piRNAs is deregulated in HNSCC, and changes with both HPV status and type. Importantly, a five-piRNA signature is able to delineate a subset of HPV-positive HNSCC patients with poor outcome, highlighting the potential utility of piRNAs in patient management.